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Abstract 

We show that as n changes, the characteristic polynomial of the n x n random 
matrix with i.i.d. complex Gaussian entries can be described recursively through a 
process analogous to Polya's urn scheme. As a result, we get a random analytic func- 
tion in the limit, which is given by a mixture of Gaussian analytic functions. This 
gives another reason why the zeros of Gaussian analytic functions and the Ginibre 
ensemble exhibit similar local repulsion, but different global behavior. Our approach 
gives new explicit formulas for the limiting analytic function. 
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Ginibre points and their random intensity measure 



1 Introduction 

In studies of random point sets in the complex plane, two canonical examples have emerged, 
sharing some features but differing in others. 

The first is the infinite Ginibre ensemble, attained as the limit of the finite Ginibre 
ensembles, i.e. the set of eigenvalues ofnxn random matrices filled with independent 
standard complex Gaussian entries. It can be thought of as the stationary distribution of a 
system of particles performing planar Brownian motions but repelling each other through 
a drift given by their inverse distance. 

The second is the zero set of the random power series ^ anz"'/ \/n\ for independent 
standard complex Gaussians a„. It is an example of a random analytic function whose 
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values are jointly centered complex Gaussian, a Gaussian analytic function or GAF, for 



short. Sodin (2000) has shown that the intensity measure (under which the measure of a 
Borel subset of C is the expected number of zeros that fall in the set) of the zeros of such 
a function in fact determines the entire distribution of zeros. The above power series, 
called planar GAF, has intensity measure given by a constant multiple of the Lebesgue 
measure. Sodin's theorem implies that this is the only Gaussian analytic function whose 
zero set has a translation-invariant distribution on C. 

Beyond translation invariance, these two processes share some local properties. For 
example, in both cases we have 

P(two points in a fixed disk of radius e) x 

where four of the six in the exponent come from the square of the area of the disk and the 
extra two comes from quadratic repulsion. Contrast this with Poisson process which only 



has e^. This behavior is shared by all non-degenerate GAFs, see the work of Nazarov and 
Sodin ( 2010[ ) on universality. 



On the global scale, however, the two models are very different. For smooth, com- 
pactly supported : C — t- M with total integral zero we have two central limit theorems 



forGinibre: E.^i^M ^ iV(0,^||V^|| 

for planar GAF: n J^z'^i^/^) =^ N{0,c\\Aip\\'' 



2) 



where the sum is over all points of the processes. These results are due to Rider and Virag 
(2007) and Sodin and Tsirelson ( 2004[ ), respectively. 



The strikingly different central limit theorems show that the global behavior of the 
two random point processes are very different. 

The goal of this paper is to prove a theorem which clarifies why this phenomenon 
happens. We study the distributional limit of the characteristic polynomial Quiz) of A„. 
It is known (see |Girko| ( [l990| ) or |Kostlan| ( |l992| )) that |(5„(0)|, the absolute value of the de- 



terminant, has the same distribution as the product of independent Gamma(i, 1) random 
variables with i = 1, . . . ,n. Since log(Gamma(z, 1)) has variance asymptotic to 1/z, we see 
that log |Qn(0)|, centered and divided by y/logn converges in distribution to a standard 
normal random variable. 

A simple consequence of this fact is that there are no constants a„, 6„ so that the ran- 
dom variable (QniO) — an)/bn converges in law to a non-degenerate limit. In contrast, we 
prove the following. Note that the convergence in distribution here is with respect to the 
topology of uniform convergence on compact subsets of the complex plane. 
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Theorem 1. There exist positive random variables An so that the normalized characteristic poly- 
nomial Qn{z)/An converges in distribution to a random analytic function Q{z). 

Moreover, there exists a random positive definite Hermitian function K{z, w) so that given K 
the function Q{z) is a Gaussian field with covariance kernel K. Further, K is analytic in z and 
anti-analytic in w, hence Q is a GAF conditional on K. 

Thus, the limit is in fact a randomized Gaussian analytic function. Theorem [l] thus 
gives a novel link between the world of random matrices and Gaussian analytic func- 
tions. In physics, certain connections between determinantal and Gaussian structures 
are referred to as super-symmetry; Theorem [l] is a very specific, mathematically precise 
instance of a connection between these two structures. 

A Gaussian analytic function (GAF) / is a complex Gaussian field on C which is almost 
surely analytic. By the theorem of Sodin ( 2000[ ) referred to earlier, the distribution of the 



zeros of a Gaussian analytic function is determined by their intensity measure 

d^i{z) = -A log K{z, z)dC{z), (1) 

TT 



where K{z,w) = E[/(z)/(w)] is the covariance kernel and C is Lebesgue measure. For- 
mula ([TJ is a special case of the well-known Kac-Rice formulas (see sec. 2.4 of Hough, 



Krishnapur, Peres and Virag ( 2009[ )). In our setting this is just an averaged version of 



Green's formula in complex analysis. 

A direct consequence is a connection between the infinite Ginibre points and Gaussian 
analytic zero points, depicted in the front-page figure. 

Corollary 2. The infinite Ginibre point process has the same distribution as Gaussian analytic 
zero points having a randomly chosen intensity measure jj,. 

Note that this does not make the infinite Ginibre point process the zero set of a true 
Gaussian analytic function, only one with a randomized covariance kernel. This random- 
ization does not change the qualitative local behavior, but changes the global one, which 
explains the phenomenon discussed above. 

We conclude our paper by computing the first moment (everywhere) and the second 
moment (at zero) of the limiting covariance kernel normalized so that these moments 
exist. 

Real Ginibre matrices. With our methods, we also prove an analogous theorem 
for characteristic polynomials of real Ginibre matrices (having i.i.d. real standard nor- 
mal entries). However there are a few modifications as the characteristic polynomials 
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are complex-valued but the values are not complex Gaussian (as complex Gaussians are 
isotropic, by definition, for us). 



Theorem 3. Let An be the real Ginibre matrix with i.i.d. N{0, 1) entries. There exist positive 
random variables An so that the normalized characteristic polynomial Qn{z)/An converges in 
distribution to a random analytic function Q{z). 

Moreover, there exist functions K{z,w), K{z,w) so that given K,K, the function Q{-) is a 
Gaussian field with E[Q{z)Q*{w)] — K{z,w) and 'E,[Q{z)Q{w)\ — K{z,w). Further, K is 
analytic in z and anti-analytic in w while K is analytic in both z and w. Therefore, Q is a random 
analytic function whose real and imaginary parts are jointly Gaussian fields. 

For later purposes, we introduce the following notation. 
Notation. We write Y ~ Normal [m( A), L(A, /i), I/(A, to mean that the real and imag- 
inary parts of Y are jointly Gaussian fields, E[F(A)] = m(A) and E[F(A)F(/x)] = L{\^) 
and E[r(A)F(/i)] = L(A,/x). 

Observe that F is a random analytic function if and only if m(A) is analytic in A, L(A, jj) 
is analytic in A and anti-analytic in /x, and L is analytic in both A and /x. Further, F is a 
complex Gaussian field if and only if L = 0. In particular, for a GAF the third argument 
is identically zero. 

A Polya's urn scheme for characteristic polynomials. We construct the random co- 
variance kernel is via a version of Polya's urn scheme. By similarity transformations, we 
first convert the matrix into a Hessenberg (lower triangular plus an extra off-diagonal) 
form. Then, we consider characteristic polynomials of the successive minors. It turns out 
that these polynomials develop in a version of Polya's urn scheme, which we recall here 
briefly. An urn containing a black and a white ball is given, and at each time an extra 
ball is added to the urn whose color is black or white with probabilities proportional to 
number of the balls of the same color already in the urn. In short, we may write 



= 0, ^2 = 1, Xk+i = Bernoulli 



k>3 



k 

For us, the essential part of Polya's urn is that the X^+i given the events up to time /c is a 
random variable whose mean is close to the average of the previous Xk. The conclusion 
is that this average converges almost surely to a random limit. 

It turns out that a similar recursion this also holds for the random constant multiples 
Xn — Xn{X) of the characteristic polynomials Qn- We shall see that 
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where 



Mfc(A,/i) = 



X, 



k+l 



Normal 



mfc(A),Mfc(A,/i),Mfc(A,/i) 



AX,(A) 



Mfe(A,/i) 



Xi(A)X*(/i) + ... + Xfc(A)X*(/i) 



bf (Xi(A)Xi(^) + . . . + Xfc(A)Xfc(/x)) 







if 13 
if /3 



1. 
2. 



with 6^7 A; 1. 

The main feature of Polya's urn is that the parameter of the Bernoulli distribution con- 
verges to a random limit almost surely, and the samples are asymptotically independent 
from this random distribution. Similarly, in our case, the variance parameter converges 
almost surely, and given the limit, the samples are asymptotically independent. In par- 
ticular, as we shall show, in the limit they behave like a Gaussian analytic function with a 
random covariance kernel. 



2 The recursion for the characteristic polynomial 



Start with the real (/3 = 1) or complex {(3 = 2) Ginibre matrix, having i.i.d real or complex 
Gaussian entries. For the purposes of this paper, a standard complex Gaussian random 
variable will mean X + iY where X, Y are i.i.d. X(0, 1/2). 



Following the randomized application of the Lanczos algorithm as pioneered by Trot- 
ter (1984]), we conjugate the matrix A by an orthogonal /unitary block matrix as follows. 
Let 



A 



ail 



\ 



All 



so that we get 



an 



( 



O 



bO 
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On 



\ 



OAO* 



One 



OiiAiiOl, 
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If Oil is chosen depending only on b so that it rotates b into the first coordinate vector, 
then we get a matrix of the form 

' AT ^X(n-l)/3 ■ • • \ 

AT AT ■■■ AT 



V 



AT AT 



AT 



where all the entries are independent, and J\f indicates that the entry has Normal (real 
or complex) distribution and Xn indicates that the entry has the distribution of the length 
of the vector with independent standard real Gaussian entries in n dimensions. This is 
because the normal vector c has a distribution that is invariant under rotation and the 
matrix An has a distribution that is invariant under conjugation by a rotation. Repeated 
application of this procedure (where the rotation matrices are block diagonal with an 
identity matrix of increasing dimension at the top) brings M to a form 



\ 



J 



We conjugate this matrix by the reverse permutation matrix and transpose it. The eigen- 
value equation of the resulting matrix reads 







7^X(n-l)/3 



X = XX. 



(2) 



Remark 4, This reduction is analogous to the tridiagonal form of the GUE matrix obtained 
by Trotter (1984]) and that has been of much use in studying the scaling limits of eigenval- 
ues in recent years (for example, see Dumitriu and Edelman| ( 2002| ), Ramirez, Rider and 
\%Igl ( |20TT] ) and |Valk6 and Virag| ( |2009l )). 



The matrix on the left of ^ is the one whose eigenvalues we will study. Indeed, let 
A/fe denote the kth row of this matrix, with the x variable removed. Then the kth. row of 
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the eigenvalue equation is 



XX k 



(3) 



where Xk = (Xi, . . . , Xk). For any A G C, let Xi{X) = 1, and define Xk+i{X) recursively as 
the solution to the equation given above for k < n. Then, A is an eigenvalue of the matrix 
if and only if it satisfies the last equation X^(A) ■ Afn = AX„(A) or equivalently, if we solve 
for k = n also and get X„+i(A) = 0. 

But these equations are consistent as n changes, so we may define the Xk{X) for all 
A; > through the infinite version of the above matrix. Clearly Xk is a polynomial of 
degree k — 1 and hence, Xk+i{X) is a random constant times the characteristic polynomial 
of the top k X k submatrix for every k. 

We introduce the random functions 



Mfc(A,/i) = Xk{X) ■ X*M = ^X,(A)X;(/x), 



Mfc(A,/^) 



X,(A)-X,(/x) = E;=i^,(A)X,(/i) 




if /3 
if /3 



1, 
2. 



where we use both z* and z to denote the complex conjugate of z. 

Let J-'k denote the a-field generated by the first k — 1 rows of the matrix together with 
Xfc/3- From ^ note that Xk+i given Tk is a Gaussian field with mean XyfjSXk/xkp covari- 
ance structure given by /SMk/xlis ^^id l^^k/xl/^- Moreover, we have 



Mk+i-Mk = Xk+i{X)Xk+i{ij), Mk+i-Mk 



Xk+l{X)Xk+l{^l) a 13 = I, 







if /3 = 2. 



Thus the evolution of Xk can be summarized as a randomized recursion with Xi = 1 and 



X, 



fc+i 



{Normal 
Normal 



AXfc Xi(A)X*(M)+...+Xfc(A)X;(M) Xi(A)Xi(^)+...+Xfc(A)Xfc(^) 



Xfc/9 ■ 



/3^fc/3 



if/3 = 2, 
if /3 = 1. 



for A; > 1. 

This is the recursion, analogous to Polya's urn discussed in the introduction. Next, we 
establish a framework for the asymptotic analysis of such recursions. 
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3 Polya's urn in Hilbert space - theorem and examples 



The main tool for the analysis of Polya's urn schemes will be the following theorem. 

Theorem 5 (Polya's urn in Hilbert space). Let Hbe a Hilbert space, and let Ri, R2, . . . Rko-i 
be deterministic elements ofH. Let be a filtration in some probability space and for each k > ko 
assume that the H-valued random variables Rk, G J^k satisfy 

Mk = + + 
k 

||E[i?fe+i|J-fc]-A4|| = {l + \\Mk\\)0{ek) + \\Rk\\0{l/k) (4) 
ElWRk+iW^lJ'k] = 0(l + ||M,_if + ||M,f) (5) 

for some positive sequence Sk such that ekk~^ is summable. 
Then converges to some limit a.s. and in L^. 

The right hand sides of Q and ^ should simply read as "small". The specific error 
terms here are taylored to the application at hand. 

Note that without the error term, Q says that the mean of Rk+i given J^^. is equal to 
Mk = {Ri + . . . + Rk) I k. This is the setting of Polya's urn scheme. 



3.1 Example: Classical Polya's urn 

Let ?^ = M and Rx = 0, i?2 = 1. Given Rx,...,Rk let 

R -\- -\- R 

Rk+i ~ Bernoulli(Mfc), with Mk = — — -. 

k 

Then Q holds without error terms and Rl^i is bounded by 1. So Mk converges to a limit 
M almost surely. 

It is well-known that M has a Beta distribution, and the limiting distribution of Rk is 
that of Bernoulli(M) (i.e., sample M and sample from Bernoulli M). In fact, for each k, the 
distribution of Rk given M is Bernoulli(M) but in our setting this is a special phenomenon 
owing to the fact that error terms vanish. 



3.2 Example: a semi-random Hessenberg matrix 

For the next example, let D c C be a closed disk. For our next example, we will consider 
the Hilbert space of 2-variable functions from D"^ C, with the following inner product 

if, 9)= / f{x,y)g{x,y)dxdy+ / f{x,x)g{x,x)dx (6) 
Jd2 Jd 
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or, in other words, consider the usual L?' inner product with respect to the sum of the 
Lebesgue measure on and Lebesgue measure on the diagonal of D^. More precisely, 
we define this inner product on smooth functions and then take the completion to get a 
Hilbert space. We will denote the corresponding norm simply by || • ||. The following 
simple Lemma is needed. 

Lemma 6. Let X ~ Normal]^, M, M] and let m : D ^ Cbe a function. Let R = (X + m) (g) 

{X + m) and R = {X + m) (g) {X + m). Then 

max{E||i?f , E||i?f } < 8(||mf + \\m^f) + 24(1 + \D\) (^||Mf + ||Mf ) . 

Here \D\ is the area of D and the norms on m, w? are L'^{D). 
Proof For h : D ^ C, the norm ^ gives 

\\h^h*\\^ = \\h^hf = \\h\\^+\\h^f (7) 

with the norms on the right hand side are for the usual L'^{D). Using this, and triangle 
inequality we get 

\\Rf = \\X + mf + \\{X + m)Y 

< 8||Xf + 8||mf + (IIX^II + ||2Xm|| + \\m^\\y 

Since ||2Xm|p < ||X^|| + ||m^||, we get the bound 

\\Rf < 8(||Xf + llX^f ) + 8(||mf + llm^f ) 

By Cauchy-Schwarz 

||Xf < ||1||||X2|| = ^/\D\ IIX^II, 

and 

EllX^f = j E\X{z)\'^dz = r]4^p (^j M^{z,z)dz + j |Mp(2, 2) ti^^ < 3 (^||Mf + ||Mf ) 

where r]A^p = E|X^| for a standard /3-Gaussian X (i.e. 3,2 for /3 = 1, 2 respectively). Putting 
all these together, we get 

E||i?f < 8(||mf + \\mY) + 24(1 + \D\) (\\Mf + . □ 
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Let bk be deterministic positive numbers such that 
One example is to take bk = \/k. 
Consider the nested matrices 



satisfy '^k ^e^ < oo. 



f M bi 
Af Af b2 



{AT 



\ 



Af J 



where the Af refer to different i.i.d. standard real (/3 = 1) or complex (/3 = 2) normal 
random variables. 

Just as in section [2], for A G C, let Xi{X) = 1, and define Xk+i{X) recursively as the 
solution of the eigenvalue equation given by the A;th row of the matrix above (for k < n). 
In other words, let A/t denote the vector formed by the first k entries of the A;th row of the 
nested matrices, write xt(A) = (Xi(A), . . . , Xk{X)), and recursively define Xk+i{X) as the 
solution to the kth row of the eigenvalue equation 

■ xt(A) + bkXk+i{X) = XXk{X). 

Let J-'k be the sigma-field generated by the first k — 1 rows of the matrix. Define 

Rk{X,fi) = (X,®X*)(A,/i) =X,(A)X*(/i), Rk{X,i2) = (X,. ® Xfc)(A, /i) = Xfc(A)X,(/i). 



Then, given J'k we have 

Xk+i ~ Normal 
Ri + . . . + Rk 



X 



k 



r-Xk, ^Mk, ^Mk 



where 



(8) 



M, 



k 



Mk 



k '^k '^k 

{Ri + ... + Rk)/k a 13 = 1, 
if /3 = 2. 



This means that conditionally on Tk the random variable Xk+i is a Gaussian field with 
the given mean function and covariance structure. For (3 = 2 it is a Gaussian analytic 
function, while for /3 = 1 it is a random anaytic function with jointly Gaussian real and 
imaginary parts. 

In order to set up a Hilbert space "H, we first fix a closed disk D c C, and consider the 
norm ^. Let Hp = Hior (3 = 2 and Hp = H x H ior P = I (the inner product onHxH 
is of course ((m, f ),(«', f')) = {u,u'){v,v')). 



11 



First consider the case (3 = 2 where we can forget Rk, Mk- Regard Rk, as random 
variables taking values in %, and as an L^(Z})-valued random variable. Then 



E[i?fc+i(A,/i) I J-fc] = ^Mfc(A,/i) + -fXfc(A)Xfc(;u) = -^Mu{\^i) + -fi?fc(A,/i). 



'^k 

Thus, condition Q follows because 



||E[i?fc+i| J-fc]-Mfe|| < i^^^||M,|| + ||i?fc||0(lA) 

= 0{ek)\\Mk\\ + \\Rk\\0{l/k) 

To check condition (|5]), we condition on (Fk and use Lemma |6] Note the conditional mean 
and variance ([T4]> of the Gaussian field X^+i. We get 



E[||i?fc+iiri^.] < ^(l + |/^|)||M,f + -(||(AX,ff +||(AX,r||) 

"k ^k 

< c\\M,r + ^^{\\xir + \\xt\\) 

= c\\Mkf + ^\\Rk\\\ 

where c depends on the sequence bk and D only, and the last equality follows from 
Writing Rk as kMk - {k - l)Mk-i we get the required upper bound c'(||Mfc|p + ||Mfc_i|p). 

Theorem |5] implies that converges (and uniformly on D) almost surely to a limit 
M. Since local convergence for analytic functions implies sup-norm convergence, the 
limit M is analytic in its two variables. Also 

Rk ^ kMk - (fc - l)Mk-i ^ Q 
k k 

and so Xk/bk 0. The conditional law of Xk+i given J'k converges to Normal [0, M, 0]. 
As each Mfc(A, /i) is analytic in the first variable and anti- analytic in the second, the same 
holds for M and it follows that Normal [0, M, 0] is in fact a Gaussian analytic function, 
conditional upon M. 

For/3 = 1, we consider the {Rk,Rk), {Mk,Mk) as elements of Hi = HxH. The estimates 
obtained above for Rk hold also for Rk with the obvious changes, and hence conditions 
^ are easily verified for {Rk,Rk)- Consequently Theorem [s] assures the existence 
of M = lim Mk and M = lim Mk and that the conditional distribution of Xk+i given Tk 
converges to Normal [0, M, M]. 
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3.3 Example: Hessenberg matrices with independent entries 

Combining the previous arguments with the central limit theorem, we can show that 
randomly scaled characteristic functions of the following Hessenberg matrices converge 
almost surely to a random analytic function which is a mixture of Gaussian analytic func- 
tions. 

( X hi \ 

X X 62 

\X ... X j 

where X are independent, identically distributed mean zero and finite fourth moment, 
and hk = \/k + 0(A;^/^~^). This condition could be significantly weakened (as long as the 
rate of growth of hk's is not too slow and not too fast); but this is not the main topic in the 
paper. 

We have universality in the sense that the limit is a randomized Gaussian analytic 
function, even when the X are not Gaussian. In particular, the local behavior of the zeros 
is universal; we expect that the probability of two zeros in an disk of radius e decays like 

e\ 

The Central Limit Theorem holds because the actual conditional covariance matrix has 
to be scaled by k to converge - in particular, a given characteristic function value has to 
use more and more of the independent X-es. Therefore it must be asymptotically normal; 
we omit the details. 

4 Polya's urn in Hilbert space - proof 

The goal of this section is to prove Theorem |5] 

We recall a few facts about Hilbert-space valued random variables and martingales. 

Fact 7. In the following, X„ will he a Hilbert-space valued random sequence. 

(i) . If Xn is a martingale, and Yl°^=i E||X„+i — X„|p is finite, then there exists a limit X so that 

||X„ — X|| — )■ almost surely and in L^{Q). 

(ii) . If E||X„ — X|p is finite, then ||X„ — X|| — )■ almost surely and in L^(n). 
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(in). If ^ (E||X„|p) is finite, then there exists a limit S so that \\S — X]fc=i ^k\\ —^0 almost 
surely and in L'^iVl). 



Proof, (ii) The claim implies that Ylm ll^n P finite almost surely (since its expected value 



is), so II X„ 1 1 -> almost surely. 



(iii) By the triangle inequality, we see that X„ is a Cauchy sequence in LriVt) and there- 



fore it converges to some S. Further, 

n oo 

\\s-Y,Xk\\ < Yl ii^'^'ii 

k=l k=n+l 

the latter being a monotone sequence that converges to in expectation, hence almost 
surely. □ 

We are now ready for the main proof. 

Proof of Theorem^ Let 

be the Doob decomposition of the process M^; the Ak are martingale increments and the 
Bk are increments of a predictable process. It suffices to show that ^ Ak and ^ Bk con- 
verge. We have 

Rk+l — E|i?fe_|_i I J^k] 



Ak 



k + 1 



which gives 

mAkinTk 



^' ^ _^ -^y (Elll-Rfe+lf I J^k] - llEli^fc+ilJ'fc]^ ) 



- ^^^Yj2E[ll^mfl'^d = ;^C'(l+||A4_if + ||Affcf) (9) 



and 



Bk = , \\Bk\\ < (1 + ||M,||)0(5.A) + \\Rk\\0{l/e) (10) 

k + 1 

First we will show that E|| M^p is bounded. Then it will follow from Fact 7|(i) that ^ Ak 



converges in L and a.s., because it is an L -bounded martingale. Then we will show that 



^(E||5fcf )i/2 is finite. Factmpi) then implies that Bi + ... + Bkis Cauchy in and 



E \\Bk\\ is finite. Thus J2 Bk converges a.s. and in L^. We write 

E||Mfc+if = E||M,,f + Epfclp + E\\Bkf + 2Re {E(iV4, Bk)} + 2E(Mfc + Bk, Ak) (11) 
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the last term vanishes, because E[Afc| J^^] = and M^, E J-'k- By (10) we have 



E\\Bkf < {l + E\\Mkf)0{el/e) + E\\RkfOil/k'') 
< (1 + E||Mfc_if + E\\Mkf)0{el/k^ + 



(12) 



The last inequality follows from the bound ([5]). Formula ( [T2| and Cauchy-Schwarz imply 
that if E|| Mfclp is bounded then ^ E||i?fc|| is finite. 
Further, two applications of Cauchy-Schwarz give 

\E{Mk,Bk)\ < E\\Mk\\\\Bk\\ 

< (E||M,fE||5,f)^/^ 

< (l + E||Mfef + E||M,_if)0(£fc/fc + l/fc2) (13) 

where the last inequality follows from ( 12 1. Using ( [12| and ( 13 1 in ( 11 1 we finally get 

E||Mfc+i||2 < E\\Mkf + (1 + E||M,f + E||M,_if )0(£fcA + l/P), 



so with the notation 



we have 



% = max EllM^ 

i<e<k 



which shows 



Vk+i < yk + il + yk + yk-i)0{ek/k + l/k'^) 
< ykil + Oisk/k + l/k^)) 

k 

yk<yil[{l + 0{ee/e+l/f)) 



in particular E||Mfc|p is bounded, completing the proof. 



□ 



5 Analysis of the Ginibre ensembles 

In this section we prove Theorem[T]and Theorem|3| By the Lanczos algorithm as explained 
in Section |2| it suffices to show that the characteristic polynomials of the Hessenberg 
matrices on the left side of ^ converge in distribution to a random analytic function 
Q, and that there exists a random positive definite Hermitian function K(z,w) so that 
given K the function Q{z) is a Gaussian field with covariance kernel K. 
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If we condition on all bl := |x^fe variables on the first super diagonal, then the matrix 



reduces to the semi-random Hessenberg matrix of Example 3.2 We leave it to the reader 



to check that k~'^\bl — k\ is summable, a.s., so that the conclusions of Example 3.2 apply. 
Thus, conditional on all bl, the limit Q of Qn (in distribution) exists, and is a mixture 
of Gaussian analytic functions. This implies that unconditionally also, Q is a mixture of 
Gaussian analytic functions. 

But we want to be able to study the properties of the random covariance kernel. Hence 
we directly analyze the Hessenberg matrix with the chi-squared variables in place, and 
get the same conclusions. 



Proof of Theorem^ Consider the nested matrices 

/ A/" 61 

Af Af b2 



{AT 



\ 



where the M refer to different real (/3 = 1) or complex {(3 = 2) normal random variables 
and b\ to Gamma (/c/3/2, 13/2) random variables, and all entries are independent. Note that 
Gamma(A;/3/2, (3/2) is the distribution of the length-squared of a fc-dimensional standard 
real or complex Gaussian vector, for (3 = 1,2 respectively. 

Define A/t as the vector formed by the first k entries of the /cth row of the nested 



matrices. As in Example 3.2 for A G C, let Xi{X) = 1, and for k < n, recursively define 
Xfc+i(A) as the solution to 

A^-xt(A) + 6fcXfc+i(A) = AX,(A) 



where 



Xfc(A) = (Xi(A),...,X,(A)). 



Then X^+i is the characteristic function of the top k k principal submatrix. 
We note that 

(3 = 2, k> 2 
/3 = 1, A; > 3 
and 

/3 = 2, A; > 3 




1 



E67 



(fc-l)(fc-2) ' 
1 

(fc-2)(fc-4) ' 
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(3 = 1, k>5 



And in light of this we will set = 1 + 4//3. Let J-'k be the sigma-field generated by the 
first k rows of the matrix. Define 



/?,(A,/i) = (X,®X,*)(A,/i) :=Xfc(A)Xfe(/x), i?,(A,/i) = {Xk®Xk){X,i2) := X,(A)X,(/i). 
Then, given and bk we have 

Xfc+i ~ Normal f^X^, f^M^, f^M^ , where (14) 

A; 

This means that conditionally on J'k and bk the random variable Xjt+i is a Gaussian field 
with the given mean and covariance. 

In order to set up a Hilbert space, we first fix a closed disk D c C Recall the Hilbert 
space v. defined in ^ and the spaces I-L2 = and T-Li = V. x V. that were introduced in 
Regard Rk, Mk, Rk, Mk as -valued random variables, and Xk as an L'^{D)- 




section 



3.2 



valued random variable. 

Theorem |5] will be applied to the 'H2-valued random variables Rk for (3 = 2 and to the 
"Hi-valued random variables {Rk, Rk) for (3 = 1. Then ^ holds as 

E[Rk+i I Tk] = nbf] {kMk + XjlXkXl) = j^Mk + ^^Rk (15) 

where A, /i are the arguments of Rk+i- Thus 

\\E[Rk+i\Fk]-Mk\\ = \\Mk\\0{llk) + \\Rk\\0{llk) 



which proves condition f or /3 = 2. Similarly, \\E[Rk+i\Fk] - Mk\\ = \\Mk\\0{l/k) + 
\\Rk\\0{l/k) from which we get 

\\E[{Rk+i,Rk+i)\J'k] - {Mk,Mk)\\ = \\{Mk,Mk)\\0{l/k) + \\{Rk,Rk)\\0{l/k) 

which proves condition Q for (3 = 1. 

We now proceed to check condition First, by Lemma |6] we have 

E[\\Rk^,r\Tk,bk] < — (l + |D|)||M,|r + -(||(AX,)2f + ||(AX,)^||) 
< ^(A:^l|M.f +||X,^f + ||X,^ 
= ^ {k'WMkr + WRkin , 
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where c depends on D only, and the last equality follows from Writing as kM^ — 
{k — l)Mk-i we get the upper bound 

E[||i?fc+if I < ^(llMfcf + ||Mfe_if). 

we complete checking condition ^ for /3 = 2 by noting that 

= 0{\\Mkf+\\Mk-iV) 

For /3 = 1, we make similar computation for Rk+i arid combine it with the above to verify 
condition 

Thus, Theorem [s] implies that (M^, Mk) converges almost surely to a limit (M, M). Of 
course M = for (3 = 2. Since local convergence for analytic functions implies sup- 
norm convergence, the limit M is analytic in the first variable and anti-analytic in the 
second variable while M is analytic in both its arguments. Also 

Rk kMk -{k- l)Mfc_i ^ 

T = k 

and so Xk/bk 0. Thus, the conditional distribution of Xk+i given converges to 
Normal[0, M, M]. For (3 = 2 this is a Gaussian analytic function with random covariance 
kernel M while for (3 = 1 it is a random analytic function with jointly Gaussian real and 
imaginary parts. □ 



6 The mean of the random covariance kernel 

We have shown that the limit of characteristic polynomials is a centered Gaussian analytic 
function with random covariance kernels M, M. Now we try to understand the distribu- 
tion of M itself. We are able to calculate the first two moments. We show the calculation 
for M only, which is sufficient for /3 = 2. 



The expectation of M(A, n) : From equation (15 1 we have 



E[Rk+i I J-fc] = E[bf] iRi + ... + Rk + XflR. 



k 



SO we would like to set to be the ERk to get the recursion 

ri + . . . + rk_+ Xjlrk 

2 
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'^fc+l = 2 (16) 



but this is invalid, because E[i?fc] = oo simply because E[6^ ^] = oo for k < 2/j3. Therefore 
we set 

rk = E[blblRk] for/3 = l, = E[6?i?fc] for /3 = 2. 



These are finite and the recursion (16) above holds for these quantities, for k > 2/(3. 



Fix z = Xfx and use ( 16 1 to write 



ri + . . . + rfc = (A; - 2/ (3)rk+i - zr^ 
ri + ... + Tk-i = {k-1 - 2/ (3)rk - zr^-i 



Taking differences we get 



z 

Tk+i -rk = 2"l^fc ~ ^^=-1 



(^3 - '"2)(fr2^ when/3 



Summing these we get 



^ " ' (r2 - ri) when (3 = 2 



= 2 : r-fc+i -ri = (r2 -ri)^— . (3 = 1: r^+i - = (rg - ra) — . (17) 



?! ^ ?! 

j=o j=o 



By direct computation we have the initial values 

1 if/3 = 2. fl + A/Z if/3 = 2. 
ri={ r2={ ra = 2 + 2A/I+ (A/Z)^ for/3 = 1. 

2 if /3 = 1. 2 + 2A7i a (3 = 1. 



Plugging these into (17 1, we get 

l + ^E,t;f + for/3 = 2. 

_2 + 2;2 + ^2 ^^^^2 ^^2 + 2z + z^e' for (3 = 1. 

as /c — 7- oo. Thus the limiting covariance kernel satisfies 



Tk+l 



E[M(A,/i)] 



1 + Xjie^f' for /3 = 2. 

2 + 2A71 + (A7Z)2e^P for /3 = 1. 



Contrast this with the planar Gaussian analytic function f{z) := En '^"7^ which has 
covariance kernel e^'^, for both real and complex i.i.d Gaussian coefficients a„ (lest this 
sound like a contradiction, one must again consider the direct second moment E[/(2)/(w)] 
which is again equal to e^*" for the real case but vanishes identically in the complex case!). 
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7 Second moment of the co variance kernel 



In this section we consider the second moment of the covariance kernel. Just like the first 
moment case, we have to multiply by random constants for this moment to exist. 

k—l 

Fix ^ = 2 and let ak = E[6f6^i?2], = B[bib^Sk], where Sk ^ RjRk- We shall get 

i=i 

recursions for these quantities, by first evaluating conditional expectations given J^^- To 
this end, set 2; = |Ap and observe that 

z'^Rl + 2k'^M^ + 4:zkMkRk 
{k-l){k-2) 

z'Rl + 2 E ■=! R'j + ^zSk + ^zRl + 4 Sj 
{k-l){k-2) 

k 

^[Sk+i\^k\ = RiY^[Rk+i\Tk\ 

i=l 

Yl k~[ ^^^^ + ^^"^ 
j=i J 

k y 

^Rj] + zSk + zR] 



k-1 

E ■=! ^1 + 2 EU + ^'^fe + ^^1 
k-1 

Multiplying by bjbl and taking expectations, we get that for k >3 



These can be rephrased as 



z{z + A)ak + 2 E ■=! + 4z/3fc + 4 E -=1 . 
(k-l)(k-2) 

Ej=i Q'j + 2 Ej=i + ^^fc + zak 



k k 
k k 

i=i j=i 



{k -l){k- 2)ak+i - z{z + 4)afc - Azjik- 
{k - l)/3jt+i - zl3k - zak- 
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For smaller k, it is possible to compute manually the following values 



k 


ak 


^k 


1 


12 





2 


6 (2 + 4^ + ^2) 


6 (1 + ^) 


3 


12 + 24 z + 24 ^2 + 8 + 


2 (6 + 9^ + 6^2 + ^3) 


4 


72 + 144^ + 144^2^96^3^33^^ + 6^5+1^6 


36 + 60^ + 48^2 + 24^3^ n 1^5 



The A; = 4 cases in fact follow from the recursions above. For A; > 4, by writing the same 
equation for /c — 1 and taking the difference we get the first two of the following equations. 
The last one follows by taking the difference of the above two. 

2ak + 4/5fe = {k-\){k- 2)afc+i -{k-2){k- 3)ak - z{z + 4)afe + ^(^ + 4)afc_i 

-4^/3fe + 4^/3fe_i. 

ak + 2/3k = {k-l)/3k+i-{k-2)/3k-z/3k + z/3k-i-z{ak-ak-i). 
2{k-l)h+i-2zPk = {k-l){k-2)ak+i- z{z + A)ak-4.zpk + '2za),. 

Using symbolic computations, we found that the coefficients of z^ is the same in all ak 
for k> j -\-2>. Similarly, the coefficients of hk — bk-i stabilize and they are exactly half the 
corresponding stable coefficients of a^. The first few coefficients are 

/ N ^,0.^0 802^2 77g^3 3799^4 9957^5 666847 11161 474659 z« 
aUz) = 72+192.+—+—+—+—+^^+^^+^^+. . . 

however, we have not been able to find a closed-form formula for the coefficients or a^{z). 
The final answer is 

lim EMfc(A,A)2 = aoo(|A|2)/2. 
Special case z = 0: The equations become 

2Q;fe + 4^fe = (A;-l)(A;-2)a;fe+i-(A;-2)(A;-3)Q;fe. 
ak + 2Pk = ik-l)Pk+,-{k-2)f5k. 

2/3k+i = {k-2)ak+i. 

Substituting the last one into the second, we get 

(k - l)(k - 2)ak+i = ak{2 + 2{k - 3) + {k - 2){k - 3)) 
which implies ak+i — ak. Then Pk — ^cts for all k >3. Thus 

k 

j=l 
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